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The recent multiple solutions problem in high energy physics has been reviewed with a more 
mathematical viewpoint. Although previously these multiple solutions are found via fit process, in 
this letter we have proved that if a sum of two coherent simple Breit-Wigner functions is used to fit 
the measured distribution, there should be two and only two non-trivial solutions, which are related 
with each other by analytical formulae. For real experimental measurements with more complex 
situations, we also provide a numerical method to extract the other solution with the already 
obtained one, and the excellent consistence between the exact solution and the fit process justifies 
this method. From our results it is clear that the physics interpretation should be very different 
depending on which solution is selected. So we suggest that all the experimental measurements with 
potential multiple solutions problem should be re-analyzed to find the other solution because the 
result is not complete if only one solution was reported. 


PACS numbers: 14.40Gx, 13.25.Gv, 13.66.Bc, 12.39.Mk 


I. INTRODUCTION 


Interference as a nature phenomenon has been ob- 
served for a very long time in situations where waves 
intersect, no matter the mediate material is water wave, 
string, sound or light. It has been studied in depth and 
also widely used in a range of physical and engineer- 
ing measurement applications, in that field Young’s slits 
interferometer, Michelson-Morley experiment and New- 
ton’s rings are famous examples. Even in some leading 
edge physics experiments, interference play a major role, 
such as optic interferometer in gravity wave detection. 
However, classic physics and quantum mechanics provide 
basically different explanations of this phenomenon. In 
classic physics, if only two meeting waves are considered 
contributing to a process, what observed is just the sum 
of amplitudes of two waves is A(x) = A(x) + A(x), 
where x are generalized coordinates could be position, 
momentum, time, energy, etc. But in quantum mechan- 
ics, their wave functions are summed to obtain the total 
amplitude (generally there are relative phases between 
them), i.e. |v(x)) = |a(x)) + |b(x)). And the experi- 
mentally measured quantities are usually proportional to 
the modulus of the amplitude squared, and thus one gen- 
erally has contribution from an interference term (a\b), 
not appears in classic physics, to the physics observable. 
Many new and fantastic features in quantum mechan- 
ics, compared with the classic one, are caused by this 
additional interference term (a|b), and the ambiguity of 
extracting information from observation is one of them. 


Usually the experimental quantities depending on |y]? 
are measured, and from which we extract the informa- 
tion of the amplitudes |a) and |b). Unlike to what in 
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classic physics, there is a square operation between the 
observable and the amplitudes, we would expect other 
solutions, |a’) and |b’), may be found in extracting am- 
plitudes from physics measurements. It’s true that the 
existing freedom on the global phase is non-relevant to 
the physics during this extraction procedure. However, 
more and more presented experimental analyses imply 
that different solutions with different relative phase an- 
gles would lead to non-trivial different physical interpre- 
tations. 


Some earlier examples reporting multiple solutions are 
in the study of the so-called Y states via initial state ra- 
diation (ISR) by the Belle experiment [1, 2]. The invari- 
ant mass distributions of rtn J/y and m+n (2S) are 
a fit with two coherent resonant terms and an incoherent 
background term. Another earlier example is the study 
of the decay dynamics of n! —> yrn mode. When the 
ntn invariant mass distribution is fitted with coherent 
sum of the p resonance and a contact term, two solutions 
are found with one solution corresponds to constructive 
interference between the two amplitudes while the other 
destructive interference [3]. Some recent examples are p- 
resented in Refs. [4, 5]. In Ref. [4] two solutions are found 
for both the branching faction measurement of ¢ > w7® 
and the p — w mixing study. In Ref. [5], four sets of so- 
lutions are found by fitting the R-values to extract the 
resonance parameters of the excited ~-family resonances, 
namely the ~(4040), (4160), and w(4415). 


However, it is notable that in Refs. [1-6] all the mul- 
tiple solutions are found via fit process. And we know 
fit method always suffers from backgrounds and limit- 
ed statistics. Then some interesting questions are raised 
naturally such as are these solutions exact solutions or 
only approximate results from the fit process? Do these 
solutions always exist or just coincidentally? How many 
are they? If one special solution has already been found, 
can the others be derived from it? Some of questions are 


explored from physics point of view in Refs. [4, 5], and 
mathematical attempts are described in Ref. [7, 8], fol- 
lowing the clue of which the present study is performed. 
Comparing with Ref. [7] in which Fourier transforma- 
tion is applied, the solution finding process is extremely 
simplified in this letter; comparing with Ref. [8], more 
general conclusions are obtained. 

In the next section, section II, at the outset, a general 
and mathematic model for the sum of two amplitudes is 
established on the basis of the known facts of physical 
analyses. If two amplitudes are both commonly used 
Breit-Wigner functions, the analytical expression for two 
solutions are obtained. Moreover, an effective approach 
is developed for acquiring the algebra equations related 
to solutions. Then many double-solutions are derived 
for the distinctive forms of amplitude functions. After 
that, we put forth the constraint relation for the ratio 
of two amplitude functions, which ensure that there will 
be double non-trivial solutions. In section III, firstly, a 
toy numerical example is used to check and confirm our 
results. Secondly, the published fit results are re-obtained 
by analytical calculation which demonstrates a numerical 
procedure to get the unknown solution from the known 
one. Finally there is a short discussion on the meaning 
about the multiple solutions problem and our suggestion. 


II. MATHEMATICAL METHODOLOGY 


If scrutinizing the relevant results of multiple solution- 
s [1-6], we can note two prominent characters: 1. all set 
of solutions have the equal goodness-of-fit; 2. although 
all parameters including the masses, the total widths, the 
partial widths and some other related parameters are al- 
lowed to float in the fit, it is observed that the only differ- 
ence between multiple solutions is the partial widths and 
the relative phase angle between amplitudes. The first 
point indicates that all solutions are mathematically e- 
quivalent while the second point implies that the main 
difference for difference solutions consists in a normal- 
ized factor and relative phase between them. In the light 
of these experimental facts, we abstract a general mathe- 
matic model for multi-solution problem. Without losing 
generality, the study that follows focuses on the case of 
two amplitude functions. 


A. Solutions for two Breit-Wigner amplitudes 


Generally, a sum of two quantum amplitudes can be 
described by a complex function e(x, 21, 22) with form 


e(z, 21, 22) = 21 g(x) + 22 f(z), (1) 
where g(x) and f(x) are both complex functions, x is a 


real variable, and z,,z2 are complex numbers. The main 
purpose of this letter is to discuss how to find non-trivial 


different series of parameters zi and z4 that satisfy 


Je(a, 21, 22)|? = jelz, 21, 29)? - (2) 


Noticed the global phase plays no role in the amplitude 
square, then we can reduce the dimension of z1 — 22 pa- 
rameter space to a z—d space in which d is real number, 
and re-write |e(x, 21, z2)|? to a more convenient form by 
defining 


le(v, 21,20)? = 5 ge) + z Fla)? 
_ a) i4 ae _ a) +2 F(a)? 
= a) E(x, z). (3) 


Here F(x) = f(x)/g(x) and E(x, z) = |1 +z F(x)|?. No- 
ticed |g(x)|? is only a multiply factor and independent of 
d and z, then it can be dropped in following discussion. 
Now we only focus on finding different series of d and z 
that keep E(x, z)/d unchanged. Denoting the real and 
imaginary parts of F(x) with R(x) and I(x), as well 
as R, and I, for z respectively. Expressing E(x, z) by 
the real and imaginary components, we obtain 


(Re + TR? + I2) —2IrI,+2RrR;+1. (4) 


Without losing generality, set d = 1 as an initial so- 
lution for convenience, so our task aims at finding all 
possible d and z’ to render E(x, z')/d = E(x, z). To spe- 
cialize our work, we consider the case when both g(x) 
and f(x) are non-relativistic Breit-Wigner (BW) ampli- 
tude functions [9] 


T 


g(x) = G- MJT, , f(z) = 


ls 
(x — My) + aT f 


where M and [I are resonance’s mass and width, re- 
spectively. This BW-form amplitude function is chosen 
because it’s universally adopted in high energy physic- 
s. With the above forms of g(x) and f(x), the real and 
imaginary components of F(x) are 


T lE Ty + (My — x)(M; = x)| 
Tar} + (My — 2)?] 


Rp = 


and 


TTM, -x)-T(Mş-— 2)] 


Ip = 
LT} + (My — «)?] 


After some algebra, we get an interesting relation: 
R24+13=aRp+blpt+c, (5) 
with 


TtT M, — M. r 
i= gts g UN é= f 
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With Eq.(5), E(x, z) is recast as 


Rr(aR?2+aI2+2R,)+Ip(bR?2 +612 -21,)+¢e(R24+-02)+1. 

(6) 
The similar expression can be obtained for E(x, 2’). No- 
tice that Rr and Ip are functions in variable space (x 
space), while R, and I, are functions in parameter space 
(z—d space), if we want for any x, E(x, 2’)/d = E(x, z), 
then the corresponding functions in parameter space as 
the coefficients of the functions in variable space should 
be equal. The requirement immediately yields: 


aR?, + al?, + 2R» 
bR2, + bI2, — 21. 
cR2, + cI2, +1 


d(aR?2 + al? +2R,), 
d(bR2 + bI? —21,) , (7) 
d(cR2 + cI? +1). 


In the light of the series of equations, it turns out that 
d must satisfy a second order equation and there’re two 
roots of it. One is the trivial solution with d = 1 and 
z’ = z correspondingly, and the other one is 


re a? + b? + 4e 
-~ (a—2R,c)? + (b+ 21c)? 


y = (ra EZD) 4 (nas BED); 


B. Some special solutions 


In this section we consider some other forms for f(x) 
and g(x). The first example is from the KLOE exper- 
iment on ete” — wr’ with both w > ata and 
w — yn°. The cross section of ete~ — wr” as a function 
of the center-of-mass energy, yz, is parameterized as 


VD onv 2 ete! O 
j Dalv) 

in Ref. [10], where onr(/z%) = oo + o'(yz — Mg) is 
the bare cross section for the non-resonant process, pa- 
rameterized as a linear function of yz; Mọ, Ty, and 
Dg = M$ — x — iMẹT¢ are the mass, the width, and the 
inverse propagator of the ¢ meson, respectively. Here z is 
a complex number which depicts the interference effect. 
Comparing with definition of E(x, z), f(x) and g(x) have 
the forms 


g(a) =—1, f(@) 


_ Mgl'¢ 
— MZ-—a2-iMglg ` 


The simple algebra yields R} + I?, = —Ip, which in turn 
gives 
R2,412,42ly = d(R2+1242L), 
Ry = dRo, (10) 
l=d. 


With the last equality d = 1, the relation R, = dR, 
implies R» = Rz, then the first equation provides the 
other non-trivial solution 


z =R,—i(I, +2). 


They are just the results acquired in Ref. [8] by another 
method. 
As the second example, we consider the from 


1 1 
which is usually used to extract the resonance informa- 
tion of w in fitting the data of ete~ > ntr. Here m 
and T indicate the mass and total decay width of the reso- 
nance. Accordingly, we obtain R} +I} = —Rpet+¢Ip(¢ = 
m/T), which in turn gives 


R477 — 2R; = d(R2+ I2 -— 2R), 
CCR, +17) - 20 = d[¢(RZ+127)—22,], (12) 
l=d. 


After some algebra, we get the other non-trivial solu- 
tion 
r 24+ (@-VR, -XT XU- Re) ~ (P - DI 
z= 
ae Se 1+0? 


(13) 
We consider a more general case, when f(x) and g(x) 
are any non-trivial functions, but their ratio must en- 
sure that the real or imaginary component of F(x) is 
constant [13]. In any of these two cases, there exist two 
solutions. Specially speaking, when F(x) = « + ih(x), 
with A(x) being a non-trivial function and « is a non- 
zero real constant, besides the trivial solution d = 1 and 
z’ = z, there exist the other solution 


1 
= 4n2(R2 + I2) + 4R: +1 
z = d[2K(R?2+/2)+ R,] + illd) (14) 


When F(x) = h(x) +i, the other non-trivial solution is 


a 1 
— 4kR(R2 + I2) — 4eI, +1 
z = R,d+ id[2«(R?2 + I?) -— L] (15) 


There are also other cases for which there exist two 
solutions, such as when Fp is a linear function of Ip, or 
vice versa, we should not discuss them in details here. 


C. Constraint for amplitude functions 


However, despite previous examples, it should be clear 
that the double-solution issue is not universal. Actually, 
it is easy to find out some forms of f(x) and g(a), in 
which no multiple-solutions can be found, g(x) = x and 
f(x) = x? + ix? is such an example. 

That’s why, although still far from the final answer, we 
want discuss what kind of constraints should be required 
in the cases when double-solutions can be found, i.e. 
what kind of amplitude functions, f(x) and g(x) in the 


preceding section, can guarantee the double-solutions. 
When we return to the study of two Breit-Wigner am- 
plitudes, we notice that the relation in Eq.(5) is crucial 
for obtaining the double-solutions. At the same time, 
this relation just provides a constraint on F(a). And it 
turns out that all the special forms with double-solutions, 
found by us, obey this requirement. Then we would like 
to extend this relation, geometrically the Argand dia- 
gram of F(x) is a circle, as a general criteria for the wave 
functions in any physics process that may take double- 
solutions. 


II. CHECK AND APPLICATION 


As across check, let’s consider an ad hoc example: the 
parameters of the two Breit-Wigner functions and one 
solution are set as: 


M, =3.0, Tg =0.4, My =2.1, Tp =0.1, z=1-%. 


Using the aforementioned method, we can find another 
solution, which is exactly repeated by fitting with mini- 
mal likely-hood method. The comparison of the results 
are shown in Table I. 


TABLE I: Comparison between exact solution and fit process 
in the case of sum of two simple Breit-Wigner functions. For 
fit, Toy MC is used to generate 10000 events data sample, the 
minimal likely-hood method is applied, all of this is realized 
in RooFit frame. 


Item | Input} Another sol. | Fit I Fit II 

d 1 0.529 = = 

R: {1 0.647 1.019 + 0.054 |0.644 + 0.040 
IL —1 1.588 —1.019 + 0.060] 1.601 + 0.028 
Ma |3.0 3.0 3.011 + 0.010 /3.011+0.010 
Tra |0.4 0.4 0.402 + 0.010 /0.402 + 0.010 
My |2.1 2.1 2.101 + 0.003 /2.101 + 0.003 
Ty, {0.1 0.1 0.101 + 0.003 /0.101 + 0.003 


This example indicates that in principle, the minimiza- 
tion can be used as a feasible approach to find the multi- 
solution from the experimental data. 

It is obvious, for two Breit-Wigner amplitudes case, 
if one solution is obtained by fitting, the other one can 
be readily and analytically obtained by applying Eq.(8). 
This definitely saves a lot of time and energy. Howev- 
er due to the complexity of the expressions in practice, 
the solution has to be obtained from following numeri- 
cal method. Firstly, we draw the F(x) in the complex 
plane to check whether it is a circle, and obtain parame- 
ters a,b,c in relation R} + Iĝ} =aRp + blr +c by three 
randomly selected points if the answer is positive. With 
a,b,c obtained numerically we can derive the other so- 
lution with Eq.(8) as before. We illustrate this method 
by the examples which are selected from the initial s- 
tate measurement at Babar and Belle [1, 12], where the 
atm w(2s) and rtr J/w invariant mass distribution- 
s are described by two coherent resonances. The cross 


sections are formulated as 
a(s) = |BW1 (5) + BW2(s) - e°}? , 


where BW, and BW, represent the two resonances and 
@ is the relative phase between them. And the Breit- 
Wigner form of a single resonance in these two papers 
is 


BW(s) = ee J/12rT e+e- B(R > fiot | PS(s) 
SEN s — M? +iMV ict PS(M) ’ 


where M is the mass of the resonance, Tot and Fete- are 
the total width and partial width to ete respectively, 
B(R — f) is the branching fraction of R decays into 
final state f, and PS(s) is the three-body decay phase 
space factor. Using the first solution as input we obtain 
the other solution as before, and the results are shown in 
TABLE II. From TABLE II it’s clear our results repeat 
what from the fit process very well, and we consider that 
is a justification of our method. 


TABLE II: Comparison of exact solutions with fit results for 
two real experimental measurements Refs.[1, 12]. 


Items BY .4,.-(R1) BTe+e- (R2) ¢ 
fit results in Ref.[1] 12.4 20.6 -111 
by our method 12.8 20.4 -111 
fit results in Ref.[12] 12.3 5.9 -74 
by our method 12.3 6.0 -74 


IV. DISCUSSION 


As been found, when the measured distribution is de- 
scribed by |g(x) + zf(x)|?/d and F(x) = f(x)/g(x) ful- 
fill the relation of Eq.(5), i.e. it’s a circle in complex 
plane, there are and only are two non-trivial solution- 
s. It’s also been proved that if f(x) and g(x) are both 
simple Breit-Wigner, this relation is exactly satisfied and 
Eq.(8) can be utilized to derive one solution from the oth- 
er obtained solution analytically. For other transmogri- 
fied Breit-Wigner functions have been considered, some 
of the forms are very complex, the relation of Eq.(5) is 
still satisfied by F(x) by numerical checks. So there will 
be double solutions for these forms also and with a, b,c 
obtained numerically the other solution can be derived 
by using the same method. The excellent consistence be- 
tween our solutions and experimental fit results justifies 
this method. 

We also notice that for both solutions, the parame- 
ters of each resonance are the same but the normaliza- 
tion factors are different. That implies the couplings to 
decay channels are different for different solutions and 
some experimental reports may not be complete if only 
one solution was reported. So we suggest any experiment 
measurement with potential multiple solutions problem 
should redo the analysis to find other solutions, and our 


method would be helpful to indicate the results. Finally, 
we should point out that from Eq.(4) we may find more 
conditions where double solutions exist, for example if 
the real or virtual component of F(x) is zero or the real 
component of F(x) is a linear function of the virtual one, 
there should be double solutions too. However, they are 
not normal in high energy physics so we don’t discuss 
it in detail here. Furthermore, only the sum of two co- 
herent amplitudes has been considered in this paper, the 
generalization to more amplitudes are still in progress. 
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